arXivil 507.08218vl [cond-mat.mes-hall] 29 Jul 2015 


Observation of droplet soliton drift resonances in a spin-transfer-torque nanocontact 

to a ferromagnetic thin film 

S. Lendmez, 1 ’^ N. Statuto, 1 D. Backes, 2 A.D. Kent, 2 and F. Macia 1 

1 Grup de Magnetisme, Departament de Fisica Fonamental, Universitat de Barcelona, Spain 
2 Department of Physics, New York University, New York, New York 10003, USA 

(Dated: July 30, 2015) 

Magnetic droplet solitons are non-linear dynamical modes that can be excited in a thin him 
with perpendicular magnetic anisotropy with a spin-transfer-torque. Although droplet solitons have 
been proved to be stable with a hysteretic response to applied currents and magnetic fields at 
low temperature, measurements at room temperature indicate less stability and reduced hysteresis 
width. Here, we report evidence of droplet soliton drift instabilities, leading to drift resonances, 
at room temperature that explains their lower stability. Micromagnetic simulations show that the 
drift instability is produced by an effective held asymmetry in the nanocontact region that can have 
different origins. 

PACS numbers: 75.78.-n, 75.78.Cd,85.75.-d,75.30.Ds 


Nanometer scale point contacts to ferromagnetic thin 
hlms with a free magnetic layer (FL) and a hxed spin- 
polarizing magnetic layer (PL)P® are known as spin 
torque nano-oscillators (STNO). In these nanocontacts, 
it is possible to excite steady state spin-precession by 
compensating the dissipation due to magnetic damp¬ 
ing with spin-transfer torques from an electrical cur¬ 
rent. A large enough current density provides suf¬ 
ficient spin-transfer-torque effect to induce magnetic 
excitation^ 5 6 . Such excitations are the building blocks 
for new applications 7 8 beyond binary computing and 
memory device^. 

In layers with perpendicular magnetic anisotropy 
(PMA) the spin-transfer torques are predicted to lead 
to dissipative droplet solitons (herafter simply refered 
as droplet solitons) 10 . The spin torque compensates 
the damping of the material, hence directly relat¬ 
ing these kind of solitons to the conservative magnon 
droplets in uniaxial (easy-axis type) ferromagnets with 
no damping 11 12 . Droplet solitons are dynam ical excita¬ 
tions that localize in the contact regiorP- 0 13 and are pre¬ 
dicted to have their magnetization almost completely re¬ 
versed relative to the him magnetization outside the con¬ 
tact [see Fig. [lj(a)]. Recent experiments have shown that 
at low temperature the predicted reversal of the mag¬ 
netization occur^Hl. However, room temperature mea¬ 
surements proved there is an abrupt threshold in both 
current and held at which excitations occur and showed 
that spin precession frequencies were below the FMR 
frequency^"^!, but there was no evidence for fully re¬ 
versed magnetization beneath the nanocontact. _ 

Here, we report observation of drift instabilities 10 ^ 
in droplet-soliton excitations at room temperature. We 
studied the effect of magnetic helds and current densi¬ 
ties on the droplet dynamics. Our modeling shows that 
drift instabilities—caused by small asymmetries in the 
system, namely a variation of either the effective held 
or the magnetic anisotropy in the nanocontact region— 
create a low-frequency dynamics (hundreds of MHz) in 
the soliton that can be detected electrically. Further, 


through simulations we have conhrmed that drift insta¬ 
bilities force the soliton to shift out of the nanocontact 
region and to annihilate while at the same time a new 
soliton emerges at the center of the contact resulting in 
a drift resonance. 
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Field 


FIG. 1. (Color Online) (a) Schematic of the STNO. An elec¬ 
trical current hows through a nanocontact (80-150 nm in di¬ 
ameter) to a thin ferromagnetic layer (free layer, FL) and 
a spin-polarizing layer (PL), (b) Schematic configuration of 
a droplet soliton centered beneath the nanocontact (left im¬ 
age) and the same soliton moved sideways (right image), (c) 
Schematic dc resistance and (d) ac response, showing the spin 
configuration of PL and FL for a nanocontact as a function 
of the out-of-plane held. 


Our samples were fabricated with a free layer (FL) 
with PMA and an in-plane hxed polarizing layer (PL). 
The FL is a 4 nm thick multilayer of CoNi, and the PL 
layer, a 10 nm thick Ni 8 oFe 2 o Permalloy (Py). A spacing 
layer of copper (Cu) magnetically decouples the two lay¬ 
ers. Si02 was used as a dielectric and Au for the contacts 
and pads. Circular contacts of different diameters were 
patterned through the dielectric by electron-beam lithog¬ 
raphy with diameters ranging from 80 to 150 nm. The 
effective anisotropy held from the CoNi multilayer was 
measured through ferromagnetic resonance spectroscopy 
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(FMR) and resulted in fi 0 H a (= hq(H k — Ms)) = 0.25 T, 
with /io the permeability of free spac£ 19 . 

The electrical response in our measurements is associ¬ 
ated with the giant magnetoresistance effect; the in-plane 
polarizer allows us to detect variations of the free layer 
magnetization in the contact region. The resistance of 
the nanocontacts depends on the relative magnetization 
orientation between the FL mpL, and PL mpL, being 
the fractional change MR = Rq(1 — mpL * mpL,)/2, with 
Ro = (Rap — Rp)/Rp, Rap,p being the resistance be¬ 
tween the device antiparallel (AP) and parallel (P) mag¬ 
netization states. An applied magnetic field perpendic¬ 
ular to the film plane tilts the Py magnetic moments 
out of the film plane, mpp • z = H/Ms for H < Ms with 
/io Ms « IT; the CoNi magnetization remains perpendic¬ 
ular to the film plane. Therefore, at large applied fields, 
H > Ms , the magnetization of the two layers forms a P 
state. 

Figure [ljc) shows schematically the point contact re¬ 
sistance as a function of the applied out-of-plane field; 
the resistance decreases linearly with the applied field 
for H < Ms and saturates when H > Ms- The mea¬ 
sured maximum overall MR of the STNO (Rq = 0.2 %) 
corresponds to twice the difference between resistance at 
zero field, where the magnetization of the PL and FL are 
orthogonal, and the resistance at a large field (H > Ms ), 
where the PL and FL are in a P state, divided by the re¬ 
sistance of the P state. The dashed red curve in Fig.[ljc) 
illustrates the expected resistance for a reversed FL mag¬ 
netization (i.e., magnetization opposite to the applied 
field). This curve is obtained by reflecting the measured 
resistance about the horizontal line, R(H = 0). 

In our experiments we detected dc variations of the 
resistance that marked the onset of droplet soliton ex¬ 
citations, the FL magnetization orientation changes at 
a threshold current while the PL magnetization remains 
fixed. We were also able to detect oscillations in the 
resistance, and thus in the voltage response, associated 
to the oscillations of the in-plane component of the FL 
magnetization. We note that the latter produces a volt¬ 
age signal for fields lower than the saturation field of the 
PL (i.e., when the polarizing layer has a component of 
magnetization in the film plane) and vanishes when the 
PL is saturated in the same direction as the FL magne¬ 
tization [see Fig. 0d)]. Finally, we also measured a low- 
frequency signal with a characteristic timescale of hun¬ 
dreds of MHz (similar to domain wall resonance^® 22 ^). 
The low-frequency signal does not vanish when the PL 
is saturated, contrary to the high-frequency precession 
of the FL magnetization [see Fig. 0d)], suggesting a 
drift resonance of the overall soliton structure beneath 
the nanocontact area. 

We first studied the onset and annihilation of droplet 
excitations through dc-measurements of resistance. We 
fixed an applied field out-of-the film plane and swept the 
applied current. These measurements show an abrupt 
increase in the resistance when the droplet forms when 
sweeping the current up — and an abrupt decrease as the 


droplet annihilates when sweeping the current down. 

Figure [2] shows the resistance curves as a function of 
the current at different fields. At each field, the cur¬ 
rent was swept up to 35 mA and then back down to 
0. Although we detected the onset of the excitations 
at fields 0.5 T < < 0.9 T, it is not until larger 

field values, /iqH > 0.9 T, that we observed hysteresis 
phenomena: while sweeping the current up, the droplet 
creation occurs at higher currents than the annihilation 
when sweeping it down. In the inset of Fig. [2] we plot a 
state map representing creation and annihilation currents 
for all measured fields. 

The hysteresis was already observed and character¬ 
ized at low temperature^. However, room temperature 
measurements showed much smaller hysteretic effectd^. 
We note here that we measured several devices and al¬ 
though the onset maps for the droplet excitations were 
almost identical, the hysteretic responses were consid¬ 
erably different: all samples showed hysteresis at larger 
fields, /io H > 0.9 T, but the size of the hysteresis varied 
between 0.5 and 5 mA. 

The size of the jumps in resistance—representing the 
difference between no excitation and droplet excitation, 
A R —are field dependent, as the angle of the PL mag¬ 
netization increases with field. The measured maximum 
change in resistance is A R = AR/Rp « 0.08 % for fields 
above the saturation of the FL magnetization (« 1 T), 
smaller than the maximum total change of Rq = 0.2 % for 
the AP configuration (being only AR/Rq « 1/3). Thus 
we conclude that the spins are not fully reversed during 
the measurement, which is a time-average measurement 
of the contact resistance. One hypothesis is that the 
magnetization of the FL precesses at an effective angle of 
about 70 degrees in the contact region; another possibil¬ 
ity is that the excitation is smaller than the nanocontact 
size, or that the excitation moves beneath th e con tact 
during the measurement time (drift instabilities 10 18 ). 
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FIG. 2. (Color Online) Measured normalized resistance R — 
R/Rp as a function of applied current for fields ranging from 
0.5 to 1.1 T. (Inset) Stability map of the droplet soliton: on 
the hysteretic area, triangles show creation of the droplet, and 
dots annihilation. 
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We have measured the high-frequency resistance signal 
of the droplet excitations. Figure [ 3 ] shows high-frequency 
spectra at a field of /iqH = 710 mT. At the onset cur¬ 
rent the characteristic frequency of the droplet excita¬ 
tion is / ~ 20.3 GHz, below the corresponding FMR fre¬ 
quency (27.5 GHz, measured in the same film), and close 
to the Zeeman frequency (y/ioTL = 19.9 GHz, 7 being the 
gyromagnetic ratio). As the current increases, the fre¬ 
quency increases (there is a blueshift) until it jumps to a 
lower frequency where it continues shifting with the same 
trend. In the inset of Fig. [3] we have plotted the power 
spectra at a fixed applied current of I = 30 mA to show 
that the quality factor of the droplet oscillations is good 
(Q ~ 2000), having peak widths of about 10 MHz. Res¬ 
onance frequencies increase with the applied field having 
a nominal value always below the FMR (between 5 and 8 
GHz). We measured similar current-dependence spectra 
at different applied fields. We note here that with in¬ 
creasing the applied field the microwave signal vanishes 
because the magnetization of the PL saturates perpen¬ 
dicular to the film plane in the same direction as the FL 
[see Fig.Qd)]. 


hysteresis in resistance. At applied fields that do not 
saturate the PL we observe different modes with increas¬ 
ing the current [see Fig. |4](a) ,(b)] and, interestingly, at 
fields that saturate the PL [see Fig. [4](c) ,(d)] the signal 
remains strong having a single peak. The magnetoresis¬ 
tance signal caused by the FL spin precession around the 
effective field vanishes when the PL is saturated; we thus 
infer that the origin of the low-frequency signal is not a 
precessing dynamics of the FL magnetization but motion 
of the whole droplet. 

Moderate in-plane fields increase the power of the 
low-frequency signal. We measured both dc and low- 
frequency MR signals as a function of the angle of the 
applied field, 0 degrees being a field perpendicular to the 
film, and obtained that the droplet annihilates and does 
not form again at angles above ±15 degrees, depending 
on the applied current. The dc resistance signal stayed at 
a higher value (corresponding to droplet state) constant 
as a function of the angle until it drops to the lower value 
(no droplet) at a certain angle; the low-frequency sig¬ 
nal increased with the angle (increased with the in plane 
field). (See Supplemental Material Sl^. 
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FIG. 3. (Color Online) High-frequency spectra as a function 
of applied current for a field of 710 mT. (Inset) Signal, in dB 
above noise, at a fixed current of I — 30 mA (white dashed 
line). The fitted data (red dashed line) shows a narrow peak 
with a FWHM = 10 MHz, and quality factor Q 2000. 

Next, we measured voltage signal at much lower fre¬ 
quencies (hundreds of MHz) in order to find oscillatory 
dynamics related to the droplet soliton as a whole ob¬ 
ject. Along with the creation of droplet excitations, we 
found a strong and broad oscillating signal at about 300 
MHz with a weak dependence on applied field and cur¬ 
rent. The low-frequency signal appears together with 
the step in resistance and we associate it with the cre¬ 
ation of the droplet excitation, see Fig. [4] All samples 
in the droplet state showed this low frequency voltage 
signal in the range of 100-800 MHz but the shape of the 
peaks were considerably different, having a well defined 
peak structure in some cases and a much broader struc¬ 
ture in others. Further, we measured hysteresis in the 
appearance of the low frequency signal along with the 


FIG. 4. (Color Online) Low frequency spectra as a function 
of applied current for fields of (a) 650 mT, (b) 889 mT, (c) 
975 mT, and (d) 1046 mT. 

We modeled the nanocontact and the droplet excita¬ 
tions with micromagnetic simulations using the open- 
source MuMax code^, and performed parallel calcula¬ 
tions with a graphics card with 2048 processing cores. 

The material parameters taken were those we deter¬ 
mined experimentally with FMR (see above). A circular 
nanocontact of 150 nm in diameter was modeled to fit 
the nominal diameter of the measured samples. We used 
a damping parameter a = 0.03 and adjusted the spin 
torque efficiency to obtain a droplet onset map similar 
to the measured in Fig. |2j We also considered the ef¬ 
fects of the Oersted fields. (Full codes are available in 
the Supplemental Material S 2 2 ^.) Additionally, we sim¬ 
ulated time frames of hundreds of nanoseconds in order 
to resolve the slow motion of the solitons. 

Our simulations show that droplet solitons form and 
annihilate at different critical current values. The hys- 
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FIG. 5. (Color Online) Time evolution of the droplet soliton 
in an applied field of 1.1 T perpendicular to the film plane 
first and with an additional in-plane field (y direction) of 0.15 
T for times t > 5 ns. The upper panels show magnetization 
maps for m Zl and m x , at particular times of the simulation. 
Images correspond to a 400 X 400 nm 2 field view. The con¬ 
tact region is outlined in black. The lower panel shows the 
time evolution of the perpendicular component of the magne¬ 
tization m z , averaged in the nanocontact area. 


teresis describes the stability of the object (with no 
temperature effects) and corresponds to a larger value 
than twice the anisotropy (2 HoH a = 0.5 T)P=! Our 
experiments showed however a much smaller hysteresis, 
«0.1T. 

Under these conditions, our simulations did not show 
drift instabilities and thus no low-frequency dynamics. 
However, we found that a small in-plane field (/jlqH ^ 0.1 
T) causes a dramatic change: droplet excitations shift in 
the direction perpendicular to the applied in-plane field 
and annihilate as they get out of the nanocontact region 
because of the damping. Immediately, a new droplet soli¬ 
ton is created beneath the nanocontactP^. This process 
has a frequency, for the simulated parameters, of about 
150 MHz and hardly depends on the out of plane field 
and applied current. 

Figure [5] shows the evolution of the droplet soliton in an 
applied field of 1.1 T perpendicular to the film plane and 
then with an additional in plane field of 0.15 T (equiva¬ 
lent to a magnetic field with an angle of 6 « 8 degrees). 
The upper panels show magnetization maps, for the com¬ 
ponents perpendicular to the plane, m z , and in the plane, 
m X: at particular times of the simulation. The lower 
panel shows the time evolution of m z in the nanocontact 
area. During the first 5 ns we apply only a perpendic¬ 
ular magnetic field of 1.1 T and observe how a droplet 
forms having all spins precessing in phase (see panels 
for t = 3 ns). At t > 5 ns we apply an in-plane field, 
0.15 T, in the y direction. This creates a drift instabil¬ 
ity, an imbalance in the precession phases that shifts the 
droplet in the x direction (perpendicular to the applied 
field) until it annihilates (see panels for t = 6 and 10 ns). 
At 12 ns it appears that the droplet has dissipated but a 


new excitation is being created (see panels for t = 12 and 
13 ns). We note that the time average of m z beneath the 
nanocontact—that is the measurable quantity using any 
dc technique—is only a 36 % of the total, or equivalent 
to a precession angle of 73 degreed. 

In order to understand why we have observed low- 
frequency dynamics even when the applied field was per¬ 
pendicular to the film we introduced asymmetric parame¬ 
ters in the simulations (see Supplemental Material Sd 23; ). 
We found that a variation in the anisotropy of only 1 % 
between the two halves of a nanocontact produces the 
same effect—with almost the same annihilation and cre¬ 
ation frequency. In general we found that any asymmetry 
in the effective field causes a drift instability and results 
in an oscillatory signal (drift resonance) of hundreds of 
MHz corresponding to annihilation and creation of the 
soliton excitation. 

Although it seems counter intuitive that hysteresis can 
exist when the droplet soliton is being created and an¬ 
nihilated, hysteresis still appears. However, the size in 
field of the hysteresis cycle reduces considerably in our 
micromagnetic simulations (from ^0.6 T to ~0.1 T). The 
reason the soliton presents certain stability even in the 
presence of drift instabilities is that when a soliton moves 
away from the contact area, the magnetization beneath 
the nanocontact is still precessing at a some finite angle. 

Drift instabilities were described as a consequence of 
magnetostatic interactions between the effective dipole 
moment of the droplet and an effective field gradient, as¬ 
sociated with the Oersted field 10 or with differences in 
the anisotropy. A small in-plane field or a small gradi¬ 
ent in the effective field creates an asymmetric landscape 
that dephases the precession of magnetization between 
edges in the soliton boundaries resulting in a magnetic 
force that acts on the soliton and shifts it. Some ex¬ 
periments have shown the presence of side bands in the 
precession frequency of the solitorPr suggesting that the 
soliton might be undergoing drift instabilities. Micro- 
magnetic studies have also predicted the side bands and 
attributed them to drift instabilities 1 -^^. 

Our direct observation of the low frequency dynam¬ 
ics proves the existence of drift instabilities and explains 
why dc measurements fail to measure full magnetization 
reversal—showing on average just a fraction of the mag¬ 
netization reversed. 

We note that the measured high-frequency dynamics, 
associated with the spin precession, has a blueshift with 
the applied current, different to predictions and micro- 
magnetic results. We cannot explain it with our model 
and we attribute such effect to the appearance of an 
effective field perpendicular to the film plane due to 
the applied current (likely an Oersted field effect from 
the leads). The overall change is always smaller than 
500 MHz that would correspond to a magnetic field of 
20 mT. 

In conclusion, we have observed and measured drift 
resonances in magnetic droplet solitons and have proved 
that droplet solitons exist and are stable at room temper- 
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ature. We suggest that the drift instability is produced 
by an effective field asymmetry in the nanocontact region 
that can have different origins. 
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Supplemental Material 


SI ADDITIONAL MEASUREMENTS 

All measured samples with the same layer stack have shown almost identical results. Here we discuss measurements 
as a function of the angle of the applied field. We first show that the results for dc resistance and frequency resistance 
signal are consistent with the results presented in the main manuscript and we discuss the small differences. The 
nanocontact sample we show in this Supplementary Material has a nominal diameter identical to the sample from 
the Main manuscript but presents a slightly larger total MR of about 0.25 % and threshold currents also slightly 
higher. We attribute this changes to small variations in the nanofabrication procedure. 

Sup. Fig. [l] shows both the dc and high-frequency resistance response measured at the same time. We see at the 
onset current, a step in dc resistance and a characteristic frequency of / ~ 20 GHz at an applied field of fi^H = 687 
mT. As the current increases, the frequency slightly blueshifts and jumps maintaining the overall value around the 
20 GHz, above the Zeeman frequency of 19 GHz 




Sup. Fig. 1, Top panel: High-frequency spectra as a function of applied current for a field of 687 mT. Bottom panel: dc MR 
as a function of current for the same sample. The two measurements were taken at the same time. 


We also measured the voltage signal at lower frequencies (hundreds of MHz). Along with the creation of droplet 
excitations we measure a strong and broad oscillating signal at about 300 MHz (see Sup. Fig.[2|. All samples showing 
droplet solitons presented the low-frequency dynamics in the range of 100-800 MHz but the shape of the peaks were 
considerably different, having a well defined peak structure in some cases and a much broader structure in others. 

Next, we measure the dependence of the droplet excitation as a function of the angle of the applied field. We first 
fixed the applied field strength and at each angle we measured the low-frequency spectra and the dc resistance. Droplet 
excitations are present only when the applied field is perpendicular to the film and we found that the maximum angle 
we can tilt the field before the excitation annihilates is about 15°. What we found is that as we tilt the angle the 
low-frequency spectra becomes stronger. Sup. Fig. [3] shows both the dc resistance and the low frequency spectra of 
the nanocontact as a function of the applied field angle. We see that the low frequency spectra is larger at angles 
between 5 and 15°. Further, we observe in the measurements shown in Sup. Fig. [ 3 ] the tilt in one direction produces 
a stronger low frequency signal. We also notice that for the stronger low-frequency signal, at positive angles between 
5 and 15° the dc MR signal also decreases a bit being consistent with the fact that the low frequency dynamics lowers 
the averaged resistance value. 
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Sup. Fig. 2. Top panel: Low frequency spectra as a function of applied current for a field of 818 mT. Bottom panel: MR as a 
function of current for the same sample. The two measurements were taken at the same time. 




Sup. Fig. 3. Top panel: Low frequency spectra as a function of the angle of the applied field for a field of 912 mT. Bottom 
panel: MR as a function of the angle of the applied field for the same sample. The two measurements were taken at the same 
time for an applied current of 40 mA 














S2 ADDITIONAL MICROMAGNETIC SIMULATIONS 


We present in here additional micromagnetic simulations that prove that spatial gradients in the effective field cause 
drift instabilities in droplet solitons. We show two representative cases, i) a variation of 1 % in the perpendicular 
applied field and ii) a variation of about 1 % in the anisotropy. 

We have divided the nanocontact in half and considered the parameters from each part slightly different. In the first 
case, we applied a perpendicular field of yoH z = 1.1 T to one half and a ~ 1 % higher in the other half (yoH z = 1.11 
T). Sup. Fig. [4] shows the evolution of the droplet soliton; first, upon an applied field of 1.1 T (3 nanoseconds), 
and then with the slight variation of 1 % in half of the nanocontact. We plotted a figure similar to Fig. 5 from the 
main manuscript showing how an imbalance is created in the soliton excitation (in the precession phases). We see 
in the upper panels of Sup. Fig. [4] how the soliton excitation shifts in the y direction until it is annihilated and a 
new excitation appears. The measured low frequency is very similar to the frequency we obtained in presence of an 
in-plane field (Fig. 5 shown in main manuscript). 
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Sup. Fig. 4. Time evolution of the droplet soliton upon an applied field of 1.1 T perpendicular to the film plane (t < 3 ns), and 
with an additional field of 0.01 T in one half of the point contact, applied at t = 3 ns. The upper panels show magnetization 
maps for m z , and m x , at particular times of the simulation. Images correspond to a 400 x 400 nm 2 field view. The contact 
region is outlined in black. The lower panel shows the time evolution of the perpendicular component of the magnetization 
m z , in the nanocontact area. 


In the second case we have considered that the anisotropies of the two parts of the nanocontact differ by 1 %. Sup. 
Fig.i shows the evolution of the soliton excitation upon an applied perpendicular field of yoH z = 1.1 T. We notice 
that in this case the dynamic annihilation and creation occurs at a lower frequency of about 50 MHz. 
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Sup. Fig. 5. Time evolution of the droplet soliton upon an applied field of 1.1 T perpendicular to the film plane and with an 
anisotropy that varies 1 % in the two halves of the nanocontact. The upper panels show magnetization maps for m z , and m x , 
at particular times of the simulation. Images correspond to a 400 x 400 nm 2 field view. The contact region is outlined in blue. 
The lower panel shows the time evolution of the perpendicular component of the magnetization m z , in the nanocontact area. 


Micromagnetics code 


//Mumax Code 

// mumaxS is a GPU-accelerated micromagnetic simulation open-source software 
// developed at the DyNaMat group of Prof. Van Waeyenberge at Ghent University. 
// The mumaxS code is written and maintained by Arne Vansteenkiste. 

//GRID 
LL := 256 
Lz := 1 

SetGridSize(LL, LL, Lz) 

SetCellSize(3e-9, 3e-9, \e-9) 

//NANOCONTACT 
diam-circ := 150e-9 
r_circ := diam-circ / 2 
A-drc := pi * pow(r-circ, 2) 

DefRegion( 1, layerf 0).interseet(circle (diam-circ))) 
tableadd(m.region (1)) 

//MATERIAL PROPERTIES 
lambda = 1 //Slonczewski parameter 

epsilonprime = 0 // Slonczewski secondairy STT term 

Msat = 500e3 //Saturation 

Kul = 200e3 //Uniaxial Anisotropy 

anis U = vector(0, 0, 1) 

Aex = 10e-12 //Exchange 

alpha = 0.03 //Damping 

//OERSTED FIELDS - 

posX := 0. 
posY := 0. 

mask := newSlice(3, LL, LL, Lz) 
current := vector(0., 0., 1.) 

for i := 0; i < LL; z++ { 
for j := 0; j < LL; j++ { 








































r ;= index2coord(i, j, 0) 
r = r. sub (vector (posX, posY, 0)) 
b := vector(0, 0, 0) 
if r.lenQ >= r-drc { 

b = r. cross (current). mul(muO / (2 * pi * r.len() * r.lenQ)) 

} 

else { 

b = r. cross (current) .mul(muO / (2 * pi * r^circ * r^circ)) 

} 

for k := 0; k < Lz; k++ { 
mask.set(0, i, j, k, b.XQ) 
mask.set(l, i, j, k, b.YQ) 
mask.set(2, i, j, k, b.Z()) 

} 

} 

} 

//END OERSTED - 

//INITIAL CONDITIONS 
angle := 85. 

my := cos(angle * pi / 180) 
mz := sin(angle * pi / 180) 

Value-Bext-Z := 1.1 //external field in z direction 
Value-Bext-X := 0.15 //external field in x direction 
angle = 90. 

px := cos(angle * pi / 180) 
pz := sin(angle * pi / 180) 

Curr := -28e-3 //current in Amps 
Pol = 0.223 

TableAddVar(Curr, ”Current”, ”A”) 

//SAVING 
tableautosave(5e-13) 
auto save (m,le-9) 


//RUNNING 

fixedlayer = vector (px,0.,pz) 
m = Uniform(0, my, mz) 


j.SetRegion(l, vector(0, 0, Curr/A_circ*Lz)) 
B_ext. RemoveExtra Terms () 

B_ext. add (mask, Curr) 

B_ext = vector(0., 0, Value-BexUz) 

run(3.0e-9) 

perpendicular field) 

B_ext = vector(Value_Bext_x, 0, Value_Bext_z) 
run(5e-8) 


//current 

// Oersted fields 
//external field 
//running 3 ns 

// new external field 
//running 50 ns 




